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In the framework of two-point sum rules of NRQCD, the masses and 
couplings of doubly heavy baryons to the corresponding quark currents are 
evaluated with account of coulomb-like corrections in the system of doubly 
heavy diquark as well as the contribution of nonperturbative terms deter- 
mined by the quark, gluon, mixed condensates and the product of gluon and 
quark condensates. The introduction of nonzero light quark mass destroys the 
factorization of baryon and diquark correlators even at the perturbative level 



I. INTRODUCTION 



>■ 

Q\ , and provides the better convergency of sum rules. We estimate the difference 

Mq - Me = 100 ± 10 MeV. The ratio of baryonic constants \Z Q \ 2 /\Z~\ 2 is 
' equal to 1.3 db 0.2 indicating the violation of SU(3) flavor symmetry for the 

doubly heavy baryons. 

o 

^ ; PACS numbers: 12.39.Hg, 14.20.Lq, 14.20.Mr 

a: 

Testing the QCD dynamics of heavy quarks in various conditions provides us with a 
qualitative and quantitative knowledge that allows us to distinguish fine complex effects 
caused by the electroweak nature of CP-violation or physics beyond the Standard Model. 
The list of hadrons containing the heavy quarks as available to the experimental observations 
and measurements, was recently enriched by a new member, the long-lived 5 c -meson, in 
addition to the heavy quarkonia bb and cc as well as the mesons and baryons with a single 
heavy quark. This success of CDF Collaboration in the first observation of i? c -meson ]]] 
was based on the progress of experimental technique in the reconstruction of rare processes 
with heavy quarks by use of vertex detectors. This experience supports a hope to observe 
other rare long-lived doubly heavy hadrons. i.e. the baryons containing two heavy quarks. 
As expected they have production rate and lifetime similar to the 5 c -meson ones. 

In the present paper we investigate the two-point sum rules @] of NRQCD ||. The light 
quark-doubly heavy diquark structure of baryon leads to the definite expressions for baryonic 



*We omit the indices QQ', when they actually do not lead to misidentificaion with usual strange 
baryons. 
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currents written in terms of nonrelativistic heavy quark fields. To relate the nonrelativstic 
heavy quark correlators to the full QCD ones we need to take into account the hard gluon 
corrections by means of solving the renormalization group equation known up to the two-loop 
accuracy. 

The convergency of sum rules results are essentially improved by account for a nonzero 
light quark mass. As was mentioned in || the sum rules stability can be achieved by de- 
stroying the baryon-diquark factorization in the correlators. The convergency was obtained 
due to taking into account the nonperturbative interactions caused by higher dimension 
operators in contrast to ref. J7|, where a signigicant instability of results were observed in 
full QCD sum rules with no product of quark and gluon condensates. We show that for the 
strange Qqq> baryons this factorization is broken already in the perturbative limit which al- 
lows us to introduce a new criterion for the determination of baryon masses since we observe 
the stability of sum rules for the masses obtained from both correlators standing in front of 
two independent Lorentz structures for the spinor field of Qqq'- 

Moreover our choice of baryonic current is convenient to take into account the a s /v 
coulomb-like corrections || inside the doubly heavy diquark. 

In Section II we describe the scheme of calculation. There we define the currents and 
represent the spectral densities in the NRQCD sum rules for various operators included 
into the consideration. The numerical estimates in comparison with the values obtained in 
potential models are given in Section III. The results are summarized in Conclusion. 

II. NRQCD SUM RULES FOR DOUBLY HEAVY BARYONS 

A. Description of the method 

In order to determine the masses and coupling constants of baryons in sum rules we 
consider the two point correlators of interpolating baryon currents. The quantum numbers 
of doubly heavy diquark in the ground states are given by its spin and parity, so that j p = 1 + 
or j p = + (if the identical heavy quarks form the diquark then the scalar state j p = + 
is forbidden). Adding the light quark to form the baryon, we obtain the pair of degenerate 
states j p = \ + and f = | + for the baryonsQ E° cc , E° bc , E° bb and H**, S*^, Ej£ with the 
vector diquark, and j p = \ + for the E bc baryons with the scalar diquark. Unlike the case of 
baryons with a single heavy quark ||, there is the only independent current component for 
each ground state. We find 

Js*, = [Q iT Cr l5 Q^q k e ijk , 

Jfcfc = lQ iT Cr 7 m Qi] . 7fn7M * eyJfc| (1) 
J5, q o q = [Q tT Cr 7 n Q J ]q%k + lriQ lT C-f m Q>] ■ l m Q%k, 



^The superscript o denotes various electric charges depending on the flavor of the light quark. 
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where J^tg satisfies the spin-3/2 condition ^ n J^*« = 0. The flavor matrix r is antisymmetric 
for H^. and symmetric for Sgg and Sqq. Here T means transposition, C is the charge 
conjugation matrix. 

The matrix structure of correlator for two baryonic currents with the spin of 1/2 has the 
form 



U(w) = i / d 4 xe tpx {0\TJ(x), J(0)|0) = fF^w) + F 2 {w) 



(2) 



where w is defined by p 2 = (M. + w) 2 , M. = mq + tjiqi + m s , ttiq^qi are the heavy quark 
masses and m s is the strange quark mass. The appropriate definitions of scalar formfactors 
for the 3/2-spin baryon are given by the following: 



IVH = i / d A xe ipx (f)\T{ J^x), J„(0)}|0) = -g^F^w) + F 2 (w)] + . 



(3) 



where we do not concern for distinct Lorentz structures. The scalar correlators F can be 
evaluated in a deep euclidean region by employing the Operator Product Expansion (OPE) 
in the framework of NRQCD, 



(4) 



where Od denotes the local operator with a given dimension d: Oq = 1, O3 = {qq), O4 = 
(-G 2 ), . . . , and the functions Cd{w) are the corresponding Wilson coefficients of OPE. For 
the contribution of quark condensate operator we explore the following OPE up to the terms 
of the fourth order in x (the derivation is presented in Appendix): 



1 



(01^(^)^(0)10) = -^<%<ss> ■ 



1 + 



X' 



(ml -2m 2 ) , x\ir 2 ^G 2 )-lm 2 s {m 2 -m 2 s )) 



16 



+ 



288 



+ im s S ab x^ n£ ■ {qq) 



1_ x^_ 
48 + 242 



3m 2 



(5) 



Note that at m s 7^ the expansion of quark condensate gives contributions in both cor- 
relators in contrast with the sum rules for Sqq/ 0, where putting m s = and neglecting 
the higher condensates, the authors found the factorization of diquark correlator in F 2 and 
full baryonic correlator in F\. This fact was the physical reason for the divergency of SR 
method. 

We write down the Wilson coefficient in front of unity and quark-gluon operators by 
making use of the dispersion relation over w, 



C,(w) = - 



71 



Pd{u)duJ 



UJ — w 



(6) 



where p denotes the imaginary part in the physical region of NRQCD. 

To relate the NRQCD correlators to the real hadrons, we use the dispersion representa- 
tion for the two point function with the physical spectral density given by the appropriate 
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resonance and continuum part. The coupling constants of baryons are defined by the fol- 
lowing expressions: 

(0\J(x)\Z(ny QQ (p)) = iZ^y QQ u{v,M m) )e^, 
(0|J m (x)|S(fi)- o (p,A))=^ H( o)*<> QM m (t;,M s(n) )e^, 

where the spinor field with the four-velocity v and mass M (the mass of baryon) satisfies 
the equation $u(v , M) = u(v , M), and u m (v , M) denotes the transversal spinor. 
Then we use the nonrelativistic expressions for the physical spectral functions 

plt^) = ^ [ \Z\ 2 5^-^ (7) 

where we have performed the substitution S(p 2 — M 2 ) — > ^5(A — w), here A is the binding 
energy of baryon and M = M. + A. The nonrelativistic dispersion relation for the hadronic 
part of sum rules has the form 



/ 



uj — w 2M A — w 



(8) 



We suppose that the continuum densities starting from the threshold u> cont , is modelled 
by the NRQCD expressions. Then, in the sum rules equalizing the correlators calculated in 
NRQCD and given by the physical states, the integration above u cont cancel each other in 
two sides of sum rules. Further, we write down the correlators in the deep underthreshold 
point of w = —M. + t with t — > 0. 

Now the sum rules in the scheme of moments with respect to t can be written down as 
follows: 



-I 

k Jo 



Wcont p 12 dw M \Z\ 



{u + M) n 2M M n ' 



(9) 



where pj contains the contributions given by various operators in OPE for the corresponding 
scalar functions Fj. Introducing the following notation for the n-th moment of two point 
correlation function: 

1 r cont p(u) du . . 

M "=-«l (fik*- (10) 

we can obtain the estimates of baryon mass M, for example, as the following: 

mm = (ii) 

and the coupling is determined by the expression 

\Z[n]\ 2 = 2 ~^M n M n+ \ (12) 
where we see the dependence of sum rule results on the scheme parameter. 
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B. Calculating the spectral densities 



In this subsection we present analytical expression for the perturbative spectral functions 
in the NRQCD approximation. The evaluation of spectral densities involves the standard 
use of Cutkosky rules fLO] , with some modification motivated by NRQCD 



1 p 2 

heavy quark : -g- -> 2m ■ 5(p - (m + — )), 

p -(m + 2m 

light quark : — - — > 2ni ■ S(p 2 — m 2 ). 

pz _ m z 

We derive the spin symmetry relations for all the spectral densities due to the fact that in the 
leading order of the heavy quark effective theory the spins of heavy quarks are decoupled, 
so 



P2,^ QQ =Sp 2Q/ o Qi =3p2,n h% , (14) 

and we have the following relation for the baryon couplings in NRQCD: 

\Z n \ 2 = 3\Z n ,\ 2 = 3\Z n *\ 2 . (15) 

Using the smallness of the strange quark mass we use the following expansions in m s for the 
perturbative spectral densities standing in front of unity operator (ttiqqi = mQmQ> / (itlq + 
ttiqi) is the reduced diquark mass, M. d i q = m Q + m Q')'- 



,3/2 



= 15015^^ + MU) + msVl ' l{u) + ^ 
where we have found 

771,0 (w) = 16u 2 (429M 3 diq + 715M 2 diq u + 403M diq u 2 + 77u 3 ), 
t] 1A (uj) = 104^(231^, + 297M 2 diq uo + \2lM diq u 2 + 15w 3 ), 
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= ( ka I , 2 (3003M dtq + 9009A*^ + 9438.M 



diq U 



+4290M 2 diq u 3 + 871M diq u 4 + 77u 5 ). (17) 

The first term of this expansion reproduces the result obtained in ||. A new feature is the 
appearance of nonzero perturbative p 2 n /<> density which is proportional to m s , 



QQ 



2^2uj(m QQI ujf' 2 m s , 
P2,n'« , M = TWTZ 3 / KA TT^rW + nv/2,i + m s T] 2t2 ), (18) 

' QQ' l\JOTC' i {JVl diq + UJ) Z 



and 
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r] 2fi = 42uj{M 2 diq + 48M diq uj + 14c; 2 ), 
V 2,i = S(35M 2 diq + 28M diq uj + 5c; 2 ), 

V2,2 = {M ^ +uj)2 (105M 3 diq + 315M 2 diq u; + 27QM diq u 2 + 77a; 3 ). (19) 

The account for the coulomb-like interaction leads to the finite renormalization of the 
diquark spectral densities before the integration over the diquark invariant mass by the 
introduction of Sommerfeld factor C, so that 

PSU* = PXuark ■ C (20) 

with 



3v, 



QQ' 



1 — exp 



27ra, x 1 1 



3vr 



(21) 



where vu denotes the relative velocity of heavy quarks inside the diquark, and we have taken 
into account the color anti-triplet structure of diquark. The relative velocity is given by the 
following expression: 



^-'/^ ^-(T-^-r (22) 

where Q 2 is the square of heavy diquark four-momentum. In NRQCD we take the limit of 
low velocities, so that denoting the diquark invariant mass squared by Q 2 = (A4 diq + e) 2 , we 
find 

3v QQI ' QQ ' 2rriQQ> ' 

at e < m Q Qi. 

The corrected spectral densities are equal to 

= m laM2M^) + + 
6ti {Jvt d i q + ui) 



where 



rj° = (2M diq + uj'fu 2 
c _ 3 (2M diq + lo )u 
{M diq + u 



^ = TTa \ v ( 12M k + 2 ±M diq u + 32M 2 diq u 2 + 20M diq u 3 + 5a; 4 ). (24) 
We see that the first term again reproduces the result of H. For the p^ / we find 
c m s m 2 QQ ,(2M diq + u)ua 

P2 = 7T-TT1 — Ta (%,() + m ^ 2 ,l + ™ s %, 2 )> ( 2 5) 

2-K{M diq + io) z 

6 



r}% = {2M d iq + 

* = ou ^ + 2^*^ + ^ (26) 

l/WeKg + UJ) 

The use of these expansions numerically leads to very small deviations from the exact in- 
tegral representations of spectral densities (about 0.5%), but they are more convenient in 
calculations. 

The contribution to the moments of the spectral densities determined by the light quark 
condensate can be calculated by the exploration of (|5|) 

M$}(n) = - ^ n+ ^ l ViM diq (n + 1) + ( " + , 3 ^ V 3 M diq (n + 3) 
n! n! 

Mfhn) = V M diq (n) - ( n+ ^ ] V 2 M diq {n + 2) + . V<tM diq (n + 4), (27) 

nl nl 

where we have introduced the coefficients of expansion in x by Vi (see (|5]) and Appendix). 
The n-th moment of two-point correlator function of diquark is denoted by Ai q (n). Then 
the diquark spectral density takes the following form: 

Pdiq = QQV , (28) 

7T 

which has to be multiplied by the Sommerfeld factor C, wherein the variable e is substituted 
by to, since in this case there is no integration over the quark-diquark invariant mass. This 
corrected density is 



liq 



(29) 



and it is independent of to. 

The corrections due to the gluon condensate are given by the density 



m Q + m Q> + HmQmQ')m n ,„ v ,c 



21 ■ 2 w ^/2Tim I Q m%{Mdiq + uj) 2 



with 



to 



t£ = 84M^ + 196M dig u + l33M dig ^ + 11 

c2 2(210.M 3 + 70M 2 diQ u + 21M diq uj 2 + 3cu 3 
Via = — 



Mdiq + UJ 

G2 _ 2(210^ + 70M|^ + 21A^ 2 + 3a; 3 ) 

where we again make the expansion in m s . In the case of nonzero quark mass we get the 
nonzero density proportional to m a , 
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^2 M = oQ 7= 2 2 ; . . ; — r (%,o + ^2,1) , (32) 

with 

' 2A M diq + uo 

For the product of condensates {gq)(^-G 2 ), wherein the gluon fields are connected to the 
heavy quarks, it is convenient to compute the contribution to the two-point correlation 
function itself. We have found 

i^-gcY .) - m QQ'( m Q + m Q' + llm Q m Q') m 
2 [LJ> 2 9 x/27rm Q m Q ,(-a;)5/ 2 ' 

and we put if G (w) = 0, since we have restricted the dimension of condensate operators, 
while the corresponding term in F\ would appear in the fifth order in expansion (gj). The 
result is given in the form, which allows the analytical continuation over uo = —M. + w. 



C. Matching with full QCD 

To connect the NRQCD sum rules to the baryonic couplings in full QCD we have to use 
the relation 

J^ = ^(« S ,/isoft,/i har d)-J™ CD , 

where the coefficient JCj(a s , // so ft, /-thard) depends on the soft normalization scale // so ft- The 
/C-factor obeys the matching condition at the hard scale //hard = M-diq an d is determined 
by the anomalous dimensions of effective baryonic currents which are independent of the 
diquark spin in the leading order. They are known up to the two loop accuracy In our 



consideration we use the one-loop accuracy, since the subleading corrections in the first a s 
order are not available yet. So, 



dln/Cj(a s ,// soft ,// hard ) ^ fa 
7 = = 2^[ — ) 7 



00 

m 

(m) 



d\n(u) ^ V47T, 

v ' m=l 

(1) = ( - 2C B {3a - 3) + 3C F (a - 2)) , (35) 



where Cf = (N% — 1)/2N C , Cb = (N c + 1)/2N C for N c = 3 and a is the gauge parameter. In 
Feynman gauge a = 1, and we get = —4. So, in the leading logarithmic approximation 
and to the one-loop accuracy we find 

7 (i) 

V ( \ ( "slXardA 2/30 , Q „x 
l^j{Oi s , // so ft, //hard) = 7 T , (36) 



a s (/t 



soft J 



where /3q = 11 — 2/3Np = 9. Further, we determine the soft normalization scale for the 
NRQCD estimates by the average momentum transfer inside the doubly heavy diquark, so 



S 



that // 2 oft = 2mQQ'Tdi q , where T^q is the kinetic energy in the system of two heavy quarks, 
which is phenomenologically independent of the heavy quark flavours and approximately 
equal to 0.2 GeV (12|. Then, the coefficients Kj are equal to 

/C Ccc ~1.95, /Cn 6c ~1.52, Ka n « 1.30, (37) 

with the characteristic uncertainty about 10% basically due to the variation of hard and 
soft scale points /ihard.soft- Note that the values of K j do not change the estimates of baryon 
masses, but they are essential in the evaluation of baryon couplings. 



III. NUMERICAL RESULTS 

Evaluating the two-point sum rules, we explore the scheme of moments. We point out 
the well-known fact that an essential part of uncertainties is caused by the variation of 
heavy quark masses. Indeed, the results of sum rules for the systems containing two heavy 
quarks strongly depend on the choice of masses, and this fact allows us to pin down the 
values of masses with a high precision up to 20 MeV [|13|], so that m;, = 4.60 ± 0.02 GeV, 



m c = 1.40 ± 0.03 GeV, which are extracted from the two-point sum rules for the families of 
T and ip. To get these values we have use the correlators evaluated up to the same accuracy 
in a s , i.e. we have put the quark loop with the appropriate Sommerfeld factor. Then, the 
stability criterion for the leptonic constant of heavy quarkonium strictly fixes the heavy 
quark masses, which are close to the results of []14| |. The same sum rules are also explored to 



estimate the couplings determining the coulomb-like interactions inside the heavy quarkonia 

a s {bb) = 0.37, a s {cb) = 0.45, a s {cc) = 0.60, (38) 

since they fix the absolute normalization of corresponding leptonic constants known exper- 
imentally. 

Since the squared size of diquark is two times larger than that of the meson the effective 
coulumb constants have to be rescaled according to the equation of evolution in QCD. We 
use the one-loop evolution equation 



a s (QQ') 

So, 



a s (QQ') 



_9_ 

Lit 



a 8 (QQ') In 2 



a s (bb) = 0.45, a s (bc) = 0.58, a s (cc) = 0.85. (39) 

The values of condensates are taken in the ranges (qq) = —(0.26 4 0.27 GeV) 3 , m 2 , = 
0.75 4 0.85 GeV 2 , (^G 2 ) = (1.7 4 1.8) • 10~ 2 GeV 4 . 

The main source of uncertainties in the ratios of the baryonic couplings is the ratio of 
the condensates of the strange quark and light quark. We use (ss) / (qq) = 0.8 ± 0.2 that 



corresponds to the variations of the sum (m u + m<i)[l GeV] = 12 4 14 MeV [15]. 

So, we have described the set of parameters entering the scheme of calculations. In Figs. 

we present the results of the two-point sum rules for the masses of E bc and VL hc (the figures 
for the other baryons are similar). For the f2;, c -baryons one can observe the stability of mass 
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with respect to the changing of the moment numbers in both correlators. We suppose it is 
connected with the destroying of diquark-f2 baryon factorization in the perturbative limit 
in contrast to the S-baryons. The stability regions for F\ and F 2 are not coincide because 
the contributions of higher dimension operators become valuable at the different numbers 
of moments. However, the quantity 1/2 (Mi [n] + M^n]) has the larger stability region, and 
we explore this fact to determine the Q baryons masses as well as that of H baryons. The 
theoretical uncertainties in the f2-baryon masses are mainly determined by the difference 
between the values of baryon masses at the regions of stability. 

Then, we investigate the difference between the masses l/2((M lj Q+M 2j n) — (M 1)S +M 2i s)) 
shown in Fig. |j. In our scheme of baryon masses determination this quantity has the 
meaning of the difference between the Q and H baryon masses. It has the large region of 
stability and is determined with a good precision. We obtain 

AM = M Qbb - M Ebb = M ncc - M Hcc = M Qbc - M= bc = 100 ± 10 MeV. 

The uncertainty in the S-baryons masses are determined through the uncertainty in the 
fi-baryons masses and that of in AM. So, for the masses we find the following results: 



M Qbc = 6.89 ±0.05 GeV, M Sbc = 6.79 ± 0.06 GeV, 

M nbb = 10.09 ±0.05 GeV, M H J = 10.00 ± 0.06 GeV, (40) 

M n = 3.65 ±0.05 GeV, M« = 3.55 ± 0.06 GeV. 

The obtained values are in agreement with the calculations in the framework of nonrela- 



tivistic potential models ||16|jl7 l, though the models based on the calculation of three body 



systems with the pair interactions |T7| give slightly higher values of masses. In || the other 
method of baryon mass determination was used, since the quantities Mi a and M 2i = sepa- 
rately have no good stability in the sum rules. So, the difference of Mi — M 2 close to zero was 
stable. The use of |(Mi + M 2 ) stability criterion results in the Sqq/ masses coinciding with 
those of |6| up to 10 MeV. Figs.|5|, |5| show the dependence of baryon couplings calculated in 
the moment scheme of NRQCD sum rules. Numerically, we find 

\Z n J 2 = (10.0 ± 1.2) -10- 3 GeV 6 , \Z E J 2 = (7.2 ± 0.8) • lO" 3 GeV 6 , 

\Z n J 2 = (15.6 ± 1.6) • 10- 3 GeV 6 , |Z H J 2 = (11.6 ± 1.0) • 10~ 3 GeV 6 , (41) 

\Z U J 2 = (6.0 ±0.8) • 10~ 2 GeV 6 , |Z H J 2 = (4.2 ± 0.6) • 10~ 2 GeV 6 . 

In Fig. we present the sum rules results for the ratio of baryonic constants 
\Zn bc \ 2 /\Zs bc \ 2 . We have also found 

\Zn b f/\Z Sb f = \Z n J 2 /\Z s J 2 = \Z Q J 2 /\Z S J 2 = 1.3 ±0.2. 

The uncertainty of this result as was mentioned above is mainly connected with the pourly 
known ratio of (ss) / (qq) = 0.8 ± 0.2. 

For the sake of comparison, we derive the relation between the baryon coupling and the 
wave function of doubly heavy baryon evaluated in the framework of potential model, where 
we have used the approximation of quark-diquark factorization. So, we find 

\Z PM \= 2^1^(0) -<M0)|, (42) 
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where ^d(O) and \I/z )S (0) denote the wave functions at the origin for the doubly heavy 
diquark and light (strange) quark-diquark systems, respectively. In the approximation used, 
the values of ^(O) were calculated in [TjJ in the potential by Buchmuller-Tye [jT8| , so that 



4tt \%(0) 
4?r |*.(0) 
4tt|* cc (0) 
4tt|^ c (0) 
4tt|* 66 (0) 



0.53 GeV 3/2 , 
0.64 GeV 3/2 , 
0.53 GeV 3/2 , 
0.73 GeV 3/2 , 
1.35 GeV 3/2 . 



In the static limit of potential models, these parameters result in the estimates 



gPM 12 


= 8.8 


10~ 3 


GeV 6 , 


\Z PM 2 


= 6.0 


10~ 3 


GeV 6 


gPM 2 


= 1.6 


10~ 2 


GeV 6 , 


\Z FM 2 


= 1.1 


10~ 2 


GeV 6 


^PM|2 


= 5.6 


10~ 2 


GeV 6 , 




= 3.9 


10~ 2 


GeV 6 



(43) 



The estimates in the potential model fl43|) are close to the values obtained in the sum rules 
of NRQCD pi|). We also see that the SU(3)-flavor splitting for the baryonic constants 
\Zn\ 2 /\Zz\ 2 is determined by the ratio |\l/ s (0)| 2 /|\I//(0)| 2 = 1.45 which is in agreement with 
the sum rules result. The values obtained in the NRQCD sum rules have to be multiplied 
by the Wilson coefficients coming from the expansion of full QCD operators in terms of 
NRQCD fields, as they have been estimated by use of corresponding anomalous dimensions. 
This procedure results in the final estimates 



\Zn cc 2 


= (38 ± 5) 


10~ 3 


GeV 6 , 


\z* I 2 

1 ^cc 1 


= (27 ±3) 


10 -3 


GeV 6 


1 7 I 2 
\^n bc \ 


= (36 ± 4) 


10~ 3 


GeV 6 , 


\z Sb f 


= (27 ±3) 


10~ 3 


GeV 6 


\Zn bb \ 2 


= (10±1) 


10~ 2 


GeV 6 , 


\z Sb f 


= (70 ± 8) 


10" 3 


GeV 6 



(44) 



IV. CONCLUSION 



In this paper the NRQCD sum rules applied to the doubly heavy baryons have been 
considered. The nonrelativistic approximation for the heavy quark fields allows us to fix 
the structure of baryonic currents (the light quark-doubly heavy diquark) and to take into 
account the coulomb-like interactions inside the doubly heavy diquark. The presence of both 
the nonzero mass of light quark and the contribution of nonperturbative terms of the quark, 
gluon, mixed condensates and the product of condensates destroys the factorization of the 
correlators. This fact provides the convergency of sum rules for each correlator and allows 
us to obtain the reliable results for the masses and baryonic constants, which agree with the 
estimates in the framework of potential models. We also have calculated the mass splitting 
of Q and H doubly heavy baryons and the ratio of baryonic constants \Z^\ 2 /\Z S \ 2 . 

The authors are grateful to prof. A.K.Likhoded for stimulating discussions. 
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V. APPENDIX 



Here the derivation of expansion (151) is briefly presented. The calculations are done in 
the technique of fixed point gauge fl9|l , so we write down the expansion of quark field: 

q{x) = q(0) + x a D a q(0) + ± x a x p D a D p q{0) + 

and in the evaluation of (0|Tg J ct (o:)g*(0)|0), where i andj are the spinor indices, a, b are the 
color indices, we have to know how to get the vacuum average of type (0|-D a .../} w g(0) g(0)|0). 
The main formulae are the followings: 
the definitions of condensates 

k a (0)g?(0)) = -^5%(gg), 

5 aa ' 

{G a al3 G a a , = -^-(g aa 'gp/3' - g a p'ga'p)(G ), 

(<1 i 9 G a^ a(T af3q}o = ml(qq), 
the commutator of covariant derivatives 

[D a ,D p ) = -igG a a/3 t a , 

and the equation of motion for the spinor field 

pq = —im q q. 

Form the last two equations we derive the so-called quadratic Dirac equation, 

D\ = -m\q + a -figG^t a q. 

Now it is an easy challenge to obtain the first term in expansion (|5|). 

Since the tensor i a ...i w is the symmetric one, we may perform the symmetrization 

D a ...D w — > {D a , Du} + , 

to find the n-th term of expansion for (q(x)q(O)), which equals 

^- x a ...x UJ {q(^) D a ...D u q(0)) = — x a ...x u) (q(Q) {L> a , ...,£>„}+ q(0)). 

Note, the tensor (q(0) {D a ...D 0J }+ q(0)) is also symmetric one. 
The second term of expansion is derived from 

({D a ,D p } + q;(0) g*(0)> = -2! V 2 ■ g a ^5 m {qq), 

and the coefficient V 2 is determined by contracting the indices a, (3 and using the quadratic 
Dirac equation, 

V 2 = (m 2 - 2mJ)/192. 
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The third term can be derived from the following structure: 

({D a , Dp, D 3 }+q p (0)q v (0)) = -3! V 3 ■ ^((t^W/^ + {lp) Pn gad + {l5) m g a p){qq)- 

Then, contracting a and (3 and using of the equation of motion, the quadratic Dirac equation 
and the commutation relation, we obtain 

V 3 = -i m q (3m 2 Q /A - mJ)/576. 

This includes the evaluation of vacuum averages 

(D 2 D a q(0) g(0)>, (D p D a Dpq(Q) q(0)) and (D a D 2 q(0) q(0)). 

Considering the structure 

({D a , Dp, D s , D $ } + q p (0)q v (0)) = -4! V 4 ■ S^S^g^ g s ^ + g aS g K + g ai . gsa){qq) 

contracted over any pair of indices, we derive 

Vi = (7c 2 (a s /n G 2 ) + 3/2 m 2 q (m 2 q - mg))/3456. 

Here we evaluated the following types of vacuum expectations: 

(D 2 D 2 q(0) q(0)), {D a DpD a Dpq{ti) q(0)), (D a D 2 D a q(0) q(0)). 

Then the OPE for the quark condensate can be expressed in terms of V% by 

(?J(ar)^(0)) = S ij (qq) ■ (V 5 pv + Vtx a ^ + V 2 5 pv x 2 + V z x al a prj x 2 + V^ m x"), 

with V = 1/12, and Vi = -im q /A8. 
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FIG. 1. The Hfe c (lower curve) and Q^c (upper curve) masses obtained in the NRQCD sum rules 
from the first correlator F\. 
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FIG. 2. The Hfe c (lower curve) and Qbc (upper curve) masses obtained in the NRQCD two point 
sum rules from the second correlator F%. 
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FIG. 3. The E;, c (lower curve) and (upper curve) masses obtained by averaging the results 
from both correlators. 
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FIG. 4. The mass difference AM = MQ bc — M= 6c obtained from the results shown in Fig 
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FIG. 5. The couplings |z£ ,2) | 2 calculated in NRQCD sum rules for the formfactors F\ and 
F2 (solid and dashed lines, correspondingly) in the scheme of moments for the spectral densities. 

\Z£ R \ 2 , GeV 6 




n 



FIG. 6. The couplings \Z r 



(1.2),2 



calculated in NRQCD sum rules for the formfactors F\ and 



F2 (solid and dashed lines, correspondingly) in the scheme of moments for the spectral densities. 
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FIG. 7. The ratio r^^ra calculated in NRQCD sum rules in the scheme of moments for the 
spectral densities at (ss)/{qq) = 0.8. 
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